Appearance of Gauge Fields and Forces beyond the adiabatic approximation 
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We investigate the origin of quantum geometric phases, gauge fields and forces beyond the adia- 
batic regime. In particular, we extend the notions of geometric magnetic and electric forces discov- 
ered in studies of the Born-Oppenheimer approximation to arbitrary quantum systems described 
by matrix valued quantum Hamiltonians. The results are illustrated by several physical relevant 
f^*) ' examples. 
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I. INTRODUCTION 

A physical system can never be considered as completely isolated from the rest of the universe. For a slow (adiabatic) 
cyclic variation of its environment, the wave function of the quantum system gets an additional geometric phase factor, 
Qh known as the Berry phase In fact, the driving environment, the 'heavy' or 'slow' system, is also subject to back 
> ■ reaction from the 'light' or 'fast' system. In the context of the Born-Oppenheimer (BO) theory of molecules, the 
. back reaction of the light system leads to the appearance of a gauge field in the effective Hamiltonian for the slow 
one (the environment) The gauge field consists of a vector and a scalar potential and turns out to depend 

a -1 

hand in the coordinate space for the evolution of trapped particles [8j and on the other hand in the momentum space 
£f~) ' for the evolution of relativistic particles (^J (see also the comment in [l^] which contains related discussion on the 

' geometric forces and fast-slow motion decoupling in physical systems), 
^vq ■ Many physical systems, as discussed in this paper, display a separation of scales in terms of slow and fast degrees 
and thus share a very similar mathematical structure with molecular systems. Actually all these systems have a 
space-time evolution governed by a multicomponent Schrodinger-like equation, whose Hamiltonian is a matrix valued 
operator. It is therefore the purpose of the present note to investigate the origin of quantum gauge fields and forces 
• • , in a more general context than the BO theory, by considering the diagonalization of an arbitrary matrix valued 
quantum Hamiltonian. To be precise, by diagonalization we mean the derivation of an effective in-band Hamiltonian 
made of block-diagonal energy subspaces. For that purpose we use the results of a method developed recently 



on a quantum geometric tensor [5j. It can both induce interference phenomena and modify the dynamics through 
geometric Lorentz and electric forces The first measurements of geometric Berry forces were done, on the one 
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This approach, based on a new differential calculus on a non-commutative space where h plays the role of running 
parameter, leads to an in-band energy operator that can be obtained systematically up to arbitrary order in H. Note 
that there exists other totally different methods of diagonalization in a formal series expansion in h which uses symbols 
of operators via Weyl calculus 12 1 13$ . Particularly important for our purpose, it has been possible, for an arbitrary 



Hamiltonian iJ(K, Q) with the canonical coordinates and momentum [Qi,Kj] — ihdij, to obtain the corresponding 
diagonal representation e (k, q) to order H 2 , in terms of non-canonical coordinates and momentum (k, q) defined later 
and commutators between gauge fields. The method is quite involved for an arbitrary Hamiltonian, but simplifies 
greatly for systems whose Hamiltonian has the simple form H = T(K) + V(Q). 

This kind of Hamiltonian that we are considering in this paper allows us to discuss how geometric gauge fields and 
geometric forces arise at order H 2 in physical situations as various as Dirac and Bloch electrons in electric fields or 
Born-Oppenheimer theory. Note that first order h corrections (semi-classical) were first treated for Bloch electrons in 



14j and for Dirac electrons in [15[|16| (see [17| for a review and also [18] which contains an overview of the first-order 



Berry-phase effects and forces in 4-D space times). 
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Our approach reveals the appearance at order h 2 of a scalar gauge potential expressed in terms of two tensors. 
One is the quantum metric tensor Q [l9| , and the other one is a new tensor generalizing an additional term found in 
12j for the Born-Oppenheimer case. Another very important consequence of the Hamiltonian diagonalization is the 
appearance of gauge invariant intraband coordinates. The advantage of using these coordinates is that the diagonal 
Hamiltonian is also gauge invariant. Moreover, these coordinates fulfill a non-commutative algebra which strongly 
affects the dynamics through a Lorentz term and the gradient of a new scalar potential, generalizing thus the dynamics 
of the Born-Oppenheimer theory. 

II. HAMILTONIAN DIAGONALIZATION. 

Consider the Schrodinger equation 

iH-U- = H\*) (1) 

where for the sake of completeness the Hamiltonian is supposed time dependent and of the form 

H = T(K,t) + V(Q,i) (2) 

where it is assumed that T(K, t) is a matrix valued operator while V(Q, t) is a scalar valued one. We aim to diagonalize 
the full differential operator D = H — Kq where Kq = ihd/dt is the conjugate operator of the time [Kq, t] = iH. 
The mathematical difficulty in performing the diagonalization of D comes from the intricate entanglement of 



noncommuting operators due to the canonical relation [Qi, Kj] — ihSij. In ll[ starting with a very general but time 
independent i?(K,Q) and by considering ft as a running parameter, we related the in-band Hamiltonian UHU + = 
e (X) and the unitary transforming matrix U (X) (where X = (Q,K)) to their classical expressions through integro- 
differential operators, i.e. £ (X) = O ^e(X)^ and U (X) = N (upC^j, where in the matrices e(X) and f/(X), the 

dynamical operators X are replaced by classical commuting variables X = (Q, K). 

The only requirement of the method is therefore the knowledge of U(K) which gives the diagonal form e(X). 
Generally, these equations do not allow to find directly e (X), U (X), however, they allow us to produce the solution 
recursively in a series expansion in h. With this assumption that both e and U can be expanded in power series of 
h, we determined in 11], the explicit diagonalization of an arbitrary Hamiltonian to order h 2 . The expression of the 



effective n-th in-band energy s n greatly simplifies for Hamiltonian given by Eq. ([2]) and the result is given by Eq. ([3]) . 
But, for this simplified problem, it is interesting to give the principal steps leading to this result. In the first step, 
assuming that the unitary matrix Uo(K,t) diagonalizing T(K,£) only, that is UqTUq = Eo(K,t), is known (eo(K,i) is 
the matrix of the eigenvalues of T(K,t)), we obtain the semiclassical expression 

U DU+ = e (K,t) + V(Q + RA,t) - K Q + HA° 

with A(K,i) = HJqV-kUq and A°(K,t) = —"iUq-^Uq two non diagonal matrices. The diagonalization at the next 
order (h 2 ) is done by an unitary transformation matrix U = Uq + HUiUq. The general method of 11 1 gives an explicit 



procedure to determine the antihermitian matrix (Ui) mn — - ^rrz^m (A)mn-^ 7 ^ which removes the off-diagonal 

elements of A and A , so that Eq. ((T|) becomes iK — ^ - = e \ ^') with \^') — U \ ^) with the diagonal energy matrix e 
whose elements e n to order h 2 are 

e„ (K, q„, t) = eo,n (K,t) + V (q„, t) + hA° + h 2 <P n (3) 
where the geometric scalar potential is 

*„(Q, K, t) = -^d i d j V+M%d i Vd j V+ «' + M;°) d % V + AC (4) 
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with two gauge invariant tensors and M£ v defined as 

G^(K,t) = lj2 * {(A% m (Ai) mn + h.c.) (5) 



2 ^— 'm^n 

and 



M%»(K,t) = W ^ )m " + h.c.) (6) 

Indices ^ corresponding to /i = 0, 1, 2, 3, such that /i = is the temporal variable and fi = i = 1, 2, 3 the spatial ones. 
The tensor GJf is known as the quantum metric tensor [5| 19 1 and M^* is a new tensor generalizing an additional 



term found in [12J for the Born-Oppenheimer theory. 

In Eq. Q the operator q = Q + HA which is non diagonal has been replaced after application of U by the 
intraband coordinate q„ = Q + ha n where a n is a gauge connection usually called the Berry connection defined as 
the projection of A = iUoVkUq on the n-th eigenstate a„(K,<) = (A) lm = i (n |Vk^)- Here \n) are the eigenstates 
of the non-diagonal part of H, i.e., T(K.t) \n) — £o.„(K,i) |n). In the same manner a° = (A°) nn = —i(n nj is 
the scalar gauge potential. Note that even though it is formally indifferent to express $„ in terms of Q or q n the 
difference being of higher order in H, the introduction of the non-canonical coordinate q ra is essential to maintain the 
gauge invariance of the Hamiltonian. The formal similarity with gauge theories is evident as we can define a Berry 
curvature the n-th eigenstate (which can be degenerate) 

leading to the following non-canonical commutation relations 

[g*,«£] =*ft a e« (8) 

As usual, the curvature leads to introduce a magnetic- type vector field 9„(K,i) whose components are defined as 
®n = E }fc®n /2 an d m the same manner from the temporal component a°(K,t) an electric-type field is defined as 

E n (K,i)=-^-V< (9) 
Note that extension of the first-order (h) Berry adiabatic formalism to the 4D space-time evolution with the "magnetic" 



and "electric" geometric fields and the corresponding Hamiltonian approach has been made in 16 1 
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For certain systems, like spinning particles in ferromagnets, the electric- type field leads to spin motive forces 
(Faraday law of spin) [20|. Here the motive force is given by 

^=-£ E - dK 4/ s e "- ds ^f 

where C is a closed curve in K space, and S the surface delimited by C (cISk being the infinitesimal surface vector 
orthogonal to S), and <f> the magnetic type flux through S. The motive force is independent of the charge of the 
particle and requires only a time dependent Berry curvature. 

The Heisenberg equations of motion q n = — ih [q„, e n ] + ^ L et K = — ih [K, e n ] to the second order in h are 



V K Eo ! „-ffl„^(Kxe„-e„xK) + / l 2 VA, K = -V q „F (10) 



The dynamics of the intraband operators leads directly to a Lorcntz-typc term. The scalar potential is a consequence 
of transitions between eigenstates \n) of T and impacts the dynamics through its gradient. Working with the non- 
canonical coordinates is a short-cut to determine the dynamics of a system prepared in a second order in h eigenstate of 
the diagonalized Hamiltonian Eq. ([3]). This state will evolve in the same energy subspace n, as there are no transitions 
between these eigenlevels as far as we can neglect higher contributions in the expansion in h. In comparison, the 
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equations of motion derived from the Hamiltonian H do not seem to include a Lorentz force,and the determination 
of the " eigendynamics" can be a very difficult to achieve. An appealing example is given in where the "exact" 
slow motion of a massive neutral particle coupled to a spin is compared with the Born-Oppenheimer theory. 

We underline that the expansion in h breaks down in regions of mode conversion where £o. m — £o.n << h or for 
large values of (m {d^n) . In a mode conversion region, one can easily generalize the diagonalization of H to a block- 
diagonalization where transitions between states inside the block are allowed (off-diagonal elements), but there are 
no transitions between different blocks. The expansion in h is consistent with the adiabatic approximation (see [l3| 
and related discussion in Indeed, the semiclassical order (h) where the potential $„ is absent corresponds to 

a situation of no-transitions between the states \n) (or between states in different blocks) whereas the expansion to 
second order in h corresponds to post-adiabatic corrections 2l( where there is no transition of higher order between 
eigenlevels and eigenstates of the diagonalized Hamiltonian Eq. ([3]). 

The virtue of our approach is its full generality which sheds a new light and provides an unified description of 
several phenomena which were considered case by case. We are now going to illustrate the results by several physical 
examples. 



III. BORN-OPPENHEIMER APPROXIMATION. 

Consider the following Hamiltonian describing a fast system in interaction with an external environment 

H=lB ij P i pi + ?-+ 9 (R,r) (II) 
2 2m 

where the fast system is described by a set of dynamical variables (r, p) (not to be confused with the non-canonical 
coordinate and momentum operators) and the slow one (the environment) by coordinates (R, P). As in Q we consider 
a general kinetic energy with B, a positive definite inverse mass tensor. Applying the previous results with the 
mapping Q — > P, K — > R (and zVk — > -«Vr) we have V^(Q) — > BijPiPj/2 and we can check that the operator 

2 

T(R) corresponds to T(R) = 2— + ^(R, r). We assume that its eigenvalues E n (R) are known (or equivalently the 
matrix Uo diagonalizing T). This eigenvalues are the energy levels of the fast system for given position R of the slow 
one, such that the matrix elements in this representation reads (m\ T(R) \n) — S m . n E n (R,). In this representation the 
Hamiltonian of the slow part is non-diagonal but we can now directly apply formula Eq. Q to obtain the following 
eigenvalues of H to order h 2 in terms of the slow variables (assuming a non-degenerate spectrum for the fast system) 



E n (p n , R) = \B ijP \jP n + hH n + E n (R) (12) 

where p n = P — ih (n | Vrr) is the gauge invariant momentum of the slow system. The scalar potential Eq. Q then 
becomes 



<!>„ (R. p) = 9njp_ B .. +M ij (R) BilBjk plpk ( 13| 



with the quantum metric tensor GJf (R) = Re^ m _^ n (djn l m ) ( TO \^j n ) an( i = Re J^m^ j^T^e^f^ • 

The term G^ (R).By is the usual part of the scalar potential discussed in several circumstances [S]j§j| 4 fl, whereas 
the term M% (R)B tl B jk P l P k was found in [12j. Here we see that the Born-Oppenheimer theory can be obtained 
straightforwardly from our Hamiltonian diagonalization to order h 2 . In the same manner from Eq.(|10p we immediately 

get the Born-Oppenheimer equations of motion p„ = — Vr-E„ — §(R X 0„ — 0„ x R) — S 2 Vr$„ with Ri = Bijp J n . 
Similar equations of motion for a classical system consisting of a classical magnetic moment interacting with an 
inhomogeneous magnetic field [6] 17:] were studied in details. It was found that the Lorentz force results from a 
slight misalignment of the magnetic moment relative to the magnetic field. This corresponds to the semi-classical 
approximation. The electric force is a time average of a strong oscillatory force induced by the precession of the 
magnetic moment. This is a kind of zitterbewegung effect. 
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IV. PARTICLE IN A LINEAR EXTERNAL POTENTIAL. 



Another interesting relevant situation concerns a particle in a linear potential exemplified here by a Bloch electron 
in a constant external electric field (see also [22j and for the first semi-classical treatment see PJ]). Consider H = 
iZo(P;R-) + ^(R) with Hq the energy of a particle in a periodic potential and </?(R) = — eE.R the external electric 
perturbation (and e < the charge). Again assuming that one knows Uq diagonalizing Hq we then have UqHqIIq = 
£o,rt (k) with £o )Tl (k) the n-th energy band and k the pseudo-momentum in the absence of the external field. Or, we 
can also write (u TO (k)| iJ (P,R) \u n (k)) = 5 m , n eo yn (k) with |zt n (k)) the periodic part of the Bloch wave function (for 
a more detailed discussion see also [23j ) . 

Then, for the determination of the full eigenvalues of H we use formula ((3]) with the mapping Q — > R, so that the 
scalar gauge potential reduces to <& n (k) = e?T^E i Ej 1 and the energy eigenvalues are 

e n = e ,n (k) - eE.r n +e 2 h 2 M^ (k) E t E, (14) 

with the covariant coordinate r„ = R + ih (u n |VkU„) and M^'(k) = R.e^ m / Wfrd^iK^ l a J M ") _ Introducing the 
"magnetic field" u n (k) = JjExVA and Xn (k) = jpE.Vk$n the equations of motion are 

r„ = V k e ,„ - |(k x n„ - O n x k) + 7l 2 x n k, k = eE 

where f2 n (k) = 0„ (k) + ftu; n (k) . This shows that $„ contributes to the Lorentz term hkx Q n known as the anomalous 
velocity which is orthogonal to the applied electric field. This anomalous velocity is at the center of many recent 



experimental and theoretical works and led to the discovery the magnetic- type monopole in solids 24|. The scalar 
potential $„ contributes also to the velocity in the direction of E, through the term fi, 2 Xnk. 



V. BEYOND THE BERRY PHASE. 



The linear potential case has another interest. It allows us to also consider the fast system to derive the Berry 
phase in a different way and to get a correction term to the phase of the wave function. 



A. General results 



Indeed, consider a time dependent Hamiltonian H{t) and introduce the differential operator D = H (t) — Pq where 
Pa = ihd/dt is the conjugate of time which is treated formally as an operator such that [Po,t] = ih. The time 
dependence is due to the time evolution of some parameters x(t) describing the environment. Assuming that we 
know e„ (t) and n(t) the instantaneous eigenvalues and eigenstates of H, i.e., H{t) \n{t)) = £„ (t) \n(t)) , then D is a 
non-diagonal matrix in this representation because of the presence of the time derivative operator Pq. To transform 
the system of differential equations (Schrodinger equation) D \^(t)) = 0, which couples all components of |\&(t)) 
into a decoupled set of differential equations, we introduce a unitary transformation |W' (t)) = U(t) ^ (t)) such that 
U(t)D (t, Pq) U + (t) = A(t, Pq) is a diagonal differential operator and A(t, Po) |\&' (t)) = 0. Therefore the time evolution 
is given Since A(t) = A(t) + P is diagonal, no time ordered product is required. 

Returning back to the initial state we have 

|* (*)) = U+(t)e^ Jo' A (*) rft C/(0) |* (0)) (15) 

A system prepared in a state |A„(0)) which is an eigenstate of D, i.e., D(0) |A„(0)) = A„(0) |A„(0)) will evolve 
with A(i) and thus stays in the instantaneous eigenstates |A n (t)) of D(t) (for simplicity we assume non degenerate 
eigenvalues). In this case the wave function becomes (t)) — e~ fo A «(*) d * |A n (t)). Since eigenstates of D instead 
of H are considered, the time evolution Eq. (fT5| is exact and thus valid for adiabatic as well as for nonadiabatic 
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evolution. Therefore for a periodic motion of period T, such that |A n (T)) = |A„(0)) (single valued eigenstates), we 
have, if |* (0)) = |A„(0)) 

with /3 n = j- J* A n (t)dt — j- e n (t) dt the exact geometric phase of the system without any approximation which 



corresponds to the Aharonov-Anandan phase |25j . Indeed, these authors extended the notion of Berry geometric 
phase for cyclic adiabatic evolutions to nonadiabatic cyclic evolutions (see also [3] and references therein). 

But, in general we need an approximation scheme for the diagonalization of D and we will use the expansion to 
order h 2 . The problem of finding A„ = A n (t) + Pq is formally equivalent to the Bloch electron example discussed 
above with K — > t, Q = R -> Pq and eE = 1. We obtain from Eq. (fT4| 



A„ = (t) - ih nj +h 2 M n (i) (17) 

x (n\m)(m n) 

with M n (t) = Rc2j m _^ n j — - — _ e 1 1 . Therefore for a periodic motion of period T, Eq. (fT6|) becomes 



|*(T))=e- < T-|*(0)) (18) 



with 



\ [ e n (t)dt + i f In n) dt +h f M n (t)dt 
"■Jo Jo * ' Jo 



The phase j n appears as an expansion in power of h. The first term is the usual dynamical phase and the second one 
the geometric Berry phase independent of h and of the velocity of parameters x(t). The additional phase h J Q M n (t) dt 
of order h is apparently non- geometric as it depends oni. It cancels in the infinitely slow x — > adiabatic regime, 
which thus coincides with the semiclassical approximation. Nevertheless, it is worth noticing that the higher order 



27 



phase corrections can also be seen as geometric in the Aharonov-Anandan meaning 27j 28| , since J Q M n (t) dt 
also be presented as a contour integral in a generalized parameter space 



Quantitatively, if the system is prepared in an eigenstate |n(0)) of -ff (0), then \^(t)) is given by Eq. (flip)) with 
U = Uq + hUiUo where U\ {t) mn — i— — ^ , so that we have the following expansion up to order h: 



|*(t)) =eT/o A "(t)d*|n(t)) + fi^ (ex/i A "^ dt A mn (t) 
- eT/.' A "^ d 'A mn (0)) \m(t))+0{h 2 ) 



i(m(t)\n(t)\ 

The magnitude of transitions is then controlled by the term A mn — £ x (t)- e (/) wm ch is neglected in the strict 



adiabatic limit (m(t) n(t)j —> 0. 

In principle deviation from adiabaticity at the semi-classical level could be measured by interferometry. Consider 
a periodic two states system, and write the initial state in the eigenbase |n(0)) = |A„(0)) + hA mn (0) |A m (0)) . Then, 
after one cycle 1^ (T)) — e~ n ™ |A„(0)) +he~ z ' lm A mn (0) |A m (0)) . For an observable O which does not commute with H 
one will find in the average (T)| O |* (T)) an interference term 2^Re (A mn (0) (A„(0)| O |A m (0)) e ' l ^~^) which 
is formally equivalent to the zitterbewegung of Dirac particles. 

Note that is normalized to unity at order h only; i.e. (^(t) \*&(t)) = 1 + 0(h 2 ) so that a normalization at 

a higher order needs an expansion of U to the same order. To second order, the diagonalizing matrix U2 (t) is given 



by the procedure described in |ll| leading to U2 (t) mn = — J2k^n m ^ (t)-£^(t))(gt(")-e (t)) so *^at ^ e wave function 
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reads 

|*(t)> = I e^/o'AnW* _ ft 2 £ e x/.' 4 "W^„ m (i)yl m „(0) + tfeT^Mt)* ([/+ (t)^ + [/ 2 ( )„„) ] |n(t)> 



-ft 2 ^ eT^WXf^fO) J |m(t)) (19) 

From this result we can compute the so called fidelity defined as K^adWI 1^(0) | 2 where \^ a d(t)) correspond to the 
wave function in the adiabatic limit. Using Eq. (|19[) we obtain 

|<n|fc>| 2 |<n|fc>| 2 



\(y ad (t)\\*(t))\ 2 = l-h 2 Y. 



J n \ (e n (t) - e k (t)f (e n (0) - e k (0)) 2 



e -i(7*-7») < n I k >< k | n > +e l ( 7fc - 7 ") I < n | fc >< fc | n > 

_|,2V ^ J — (201 

fe ^ ( £ „(*)- £fc (*))( £ „(0)- Efc (0)) 1 J 

with e -*(7»-7.) = e -if[^ + ^-<t>n) + n(M k -M n) ] dt w . th ^ = ^ 

B. Example 

As a physical illustration let consider the Hamiltonian 

H = B(t)-a 

corresponding to the paradigmatic example of a spin coupled to a magnetic field. The eigenvalues of D = H (t)—ihd/dt 
at order h 2 are given by Eq. (J 171) 

A, = £l (t) + Haf+h 2 M l (t) (21) 

The Hamiltonian being diagonalized by the matrix C/q = (™' cr+a ^ w jth n = B/B, we have the two level eigenvalues 
ei (t) = B(t) (cr z ) u with 1 = 1,2, and the matrix A = —iMJo^Uq is thus given 

nn + k) xnj .cr S 



.4 U 



2(l + n z ) 



/ ri n /i ( (n+k) X n J -k 

As a consequence the scalar gauge potential a] = — i < I | I >= —ih [Uo-^Uq ) u reads a" = 2(i+n ) — and 

Mi (t) = A/; 00 = i Em/if 1 !^^ + ^-c-) can therefore be written 

2(l + n z )n — n z — ((nxn) • k) 

Mi = 5-^ '- '— (22) 

8(l + n z fei(t) 
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n x n ] .k 



The Berry phase can thus be cast in the form 

ha = -(*.)„ f (23) 

It is interesting to introduce the Euler angles n = ( sin# cos tp, sin# simp, cos 0), so that the Berry phase reads <Pi.b — 
(<J z )u 5 / £ (1 — cos#) dtp and the correction of higher order to the phase 7; is thus 

fT 2 .2 

M ; (t) eft = — / (sin 2 0y> + 6> )eft. (24) 



We retrieve the result of |27| for y> = 0. Note that, it was shown in |27j, this higher order contribution to the geometric 
phase is a part of the non-adiabatic Aharonov-Anandan geometric phase. 

Let compute the fidelity in the simple case 9 = const and tp = cut. A direct application of Eq. (|20|) gives 

, , , , , * , ,9 tjJ 2 Sill O n , . „, 

\(* ad (t)\ |*(t))| 2 = 1 r-5— sin 2 (wo - w) t + ... (25) 

■ 2 a 2 

where we introduced the notation B = hajQ. The usual criteria for adiabaticity is given by the condition s '°Ja << 1- 
At the resonance — = 1 fidelity is always one at this order of the expansion. 



VI. DIRAC PARTICLE IN AN EXTERNAL POTENTIAL. 

We will now show that our formalism can also be used for relativistic Dirac particles, which are usually treated 
with the Foldy Wouthuysen approach 29]. The semiclassical treatment of the problem was first considered in [l5| 
and 16 1 . The direct second order in H diagonalization being done in [ill ] , we can now fully appreciate how simple and 
straightforward is the application of the formalism presented here. 

The Hamiltonian is (with c = 1) 

H = a.p + /3m + V (R) (26) 

where a and /3 are the usual (4 x 4) Dirac matrices and V (R) is the external potential. The matrix diagonalizing 
the free part of the Hamiltonian Uq (a.p + j3m) Uq = (3E with E = \Jp 2 + m 2 is the usual Foldy Wouthuysen 
unitary transformation Uq = pfj^E^fypTa ■ For the Dirac particles we have two energy subspaces e± of dimension 2 
corresponding to the positive and negative energy. Now with the correspondence Q > R, K > p and formula Eq. 
© , one easily sees that the diagonal matrix can be written 

s(p,r) = /3E(p) + V(r)+l3<f> (27) 

The position operator is given by the (4 x 4) matrix r = R+ 2 E(E+m) w ^ n ^ = 1 ® cr where a are the Pauli matrices. 
The band index of the scalar potential has been transferred to the matrix /3, and we have G y = j^g 1 -* and M y = 
gjp-g 1 - 5 with the notation g 1 ^ = 5^ — so that finally we can write 



Note the presence of a new term not presented in other considerations of the Dirac equation (see for instance [30|), 
which is nonlinear in scalar potential and which stems from the new tensor M lJ . If for central potential one can 
neglect this new contribution ^diVdjV, this is not always true and for some potentials both terms in Eq. (|28|l can 
be of the same magnitude. In fact for constant electric field V = — e£.R, the first term vanishes and $ = \§r9 lJ 

In the non-relativistic limit p << m, $ becomes $ ~ (AT^+^(VT^) 2 )) + 0(h 2 p 2 /m 4 ) which gives two contri- 
butions. The first one is the usual Darwin term AV traditionally obtained as the result of the Foldy Wouthuysen 
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transformation expanded in power of 1/m. The second term (W) 2 of higher order in 1/m is usually not consid- 
ered in the Foldy Wouthuysen approach. It is also interesting to note that the external potential in the non relativistic 
limit can be expanded as V (r) sa V (R) + ^S. (W x p) + 0(H 2 p 2 /m A ) where ^jE. (W x p) is the spin-orbit 
coupling term. Therefore the Hamiltonian can be approximated as 

/ p2 p4 \ t 

\ 2m omr J 4m z 

+ -^/3 f AV+-(VF) 2 ) (29) 
8m/ \ m / 

A Born-Oppenheimer treatment of the Dirac equation where the spin is the fast variable and the momentum the 
slow one has led to the same Hamiltonian Eq. (|2T)1) but without the scalar potential [311 ]. This corresponds to 
the semiclassical approximation. The additional electric-type potential $ is a consequence of transitions between 
energy levels. This is in agreement with the usual interpretation of the physical origin of the Darwin term, the 
zitterbewegung phenomenon, whereby the electron does not move smoothly but instead undergoes extremely rapid 
small-scale fluctuations due to an interference between positive and negative energy states. 



VII. CONCLUSION. 

In this paper we investigated the origin of quantum geometric phases, gauge fields and forces beyond the adiabatic 
approximation for physical systems displaying a separation of scales in terms of slow and fast degrees. In particular 
we extended the notions of geometric magnetic and electric forces discovered in studies of the Born-Oppenheimer 
approximation. Our approach is very general and results found here have been straightforwardly applied to several 
physical systems. 
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